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ABSTRACT 



The present work investigates the parametric instability of parallel propagat- 
ing circularly polarized Alfven (pump) waves in a weakly ionized molecular cloud. 
It is shown that the relative drift between the plasma particles gives rise to the 
Hall effect resulting in the modified pump wave characteristics. Although the lin- 
earized fluid equations with periodic coefficients are difficult to solve analytically, 
it is shown that a linear transformation can remove the periodic dependence. The 
resulting linearized equations with constant coefficients are used to derive an al- 
gebraic dispersion relation. The growth rate of the parametric instability is a 
sensitive function of the amplitude of the pump wave as well as to the ratio of 
the pump and the modified dust-cyclotron frequencies. The instability is insen- 
sitive to the plasma (3. The results are applied to the molecular clouds. 

Subject headings: Waves-ISM:clouds-MHD-stars:formation 
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1. INTRODUCTION 

Molecular clouds are thought to be the nursery of star formation (Lada 1985). These 
clouds generally consist of very small (~ 10~ 4 — 1CT 8 ) fraction of ionized matter and 
yet, depending upon the level of plasma-neutral coupling, the magnetic field can play an 
important role in the dynamics of such a cloud. Thus the question of how well the magnetic 
field is coupled to the neutral matter becomes crucial for various models of cloud support 
and star formation. The collisonal momentum exchange and the ratio of the gyro to various 
collision frequencies determine this coupling. The relative drift between the plasma particles 
introduces another interesting feature to this complex scenario. For example, when ions 
and neutrals are strongly coupled, the relative drift of the "frozen-in" ions against the sea 
of neutrals may cause the diffusion of the magnetic flux reducing the flux to mass ratio in 
the star forming regions (Mestel & Spitzer 1956; Mouschovias et al. 2006). When the ions 
and neutrals are moving together and there is relative drift between the electrons and ions, 
the redistribution of energy and angular momentum in the medium may take place (Wardle 
1999). Thus, in molecular clouds, investigation of the non-ideal magnetohydrodynamic 
(MHD) effects is important. There are three regimes in which non-ideal MHD effects 
may operate, (I) ambipolar - when the magnetic field can be regarded as frozen in the 
plasma and drifts with it through the neutrals (applicable to the relatively low density, 
high ionization fraction regions), (II) Ohmic - when neutrals stops the ionized particle from 
drifting with the field (applicable to the high density, low ionization fraction regions) and 
(III) Hall when electrons are well coupled or partially coupled to the field and ions and 
grains are partially or completely decoupled from the field. 

The magnetic field can provide either static or dynamic (by magnetohydrodynamic 
waves) support to the molecular clouds (Dewar 1970; Mouschovias & Spitzer 1976; Shu et 
al. 1987; Martin et al. 1997; Fiege & Pudritz 2000). The Alfven waves are not easily 
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dissipated (Zweibel & Josafatsson 1983) making them an attractive candidate for the 
dynamic support of the cloud. The origin of the non-thermal line broadening in molecular 
cloud is attributed to the Alfvenic turbulence owing to the non-steepening of the large 
amplitude Alfven waves (Falle & Hartquist 2002; Folini et al. 2004; Mouschovias et al. 
2006). The supernova shocks, gravitational collapse or the large scale motion of the cloud 
could act as the source of such waves in the medium. 

In dense molecular cloud cores charged grains are more numerous in number than the 
electrons and ions (Wardle & Ng 1999). The ionization fraction of the plasma is strongly 
affected by the abundance and size distribution of the grains through the recombination 
process on the grain surface. In molecular clouds with densities n n > 10 3 cm~ 3 , inelastic 
collision of various charged species with the grain is important (Ciolek & Mouschovias 

1993) . At the densities relevant to cloud cores and at higher densities occurring during the 
formation of protostar and protostellar discs (densities > 10 7 — 10 11 cm -3 ), the role of the 
grains can not be neglected in the dynamics (Wardle 1999). The presence of the charged 
grains in the planetary and interstellar medium causes the excitation of very low frequency 
Alfven waves (Pilipp et al. 1987; Cramer & Vladimirov 1994; de Juli et al. 2005). These 
waves provide important physical mechanism for the transport of angular momentum and 
energy in a differentially rotating dusty disks (Wardle 1999; Pandey & Wardle 2006). 

The parametric instability of a finite amplitude circularly polarized Alfven wave 
have been studied for last several decades (Derby 1978; Goldstein 1978; Hollweg 

1994) . Recently, this problem has been studied in the single fluid MHD framework for 
a self-gravitating molecular cloud (Fukuda & Hanawa 1999). We shall leave the Jeans 
instability out of the present investigation and revisit the problem of parametric instability 
in the molecular clouds by including the charged grains and neutrals in the formulation. 
The investigation of parametric instability in a dusty medium is relatively recent phenomena 
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(Vladimirov et al. 2005; Hertzberg et al. 2003). The charged grains are assumed immobile 
and the pump wave is assumed magnetosonic or Alfvenic. A recent generalization of the 
parametric investigation to the multi-component plasma (Hertzberg et al. 2004) removes 
the limitation of the immobile dust grains. However, the dynamics of the neutrals has been 
left out in making it inapplicable to the interstellar medium. 

It is known that collisonal processes in weakly ionized plasma can strongly influence 
the nonlinear ambipolar diffusion and the ensuing current sheet formation (Cramer & 
Vladimirov 1994). The role of the plasma-dust collision may as well be important for the 
parametric instability. For example, collisions between the plasma particles and grains are 
responsible for some of the novel features in dusty plasma. The new collective behavior 
is known to exist in such plasma due to the charge fluctuations - an offshoot of collision 
(Vladimirov et al. 2005; Bhatt & Pandey 1994). Therefore, it is desirable to investigate 
the propagation of the large amplitude Alfvenwave in a weakly ionized collisonal plasma. 
We note that in dark molecular clouds, the momentum of the bulk fluid is carried by the 
neutrals and the current is carried by the ionized component. As we shall see this will 
provide some economy of description to the present problem. 

One of the difficulties in the investigation of parametric instability of circularly polarized 
Alfven wave is associated with the appearance of periodic coefficients in the linearized MHD 
equations. This restricts the normal mode analysis of the spatial dependence and thus, 
the resulting Floquet equation can only be solved in the long wavelength limit. However, 
a recently proposed linear transformation (Ruderman & Simpson 2003) reduces the 
periodic coefficients in the linearized MHD equations to the constant coefficient. Although 
the transformation by Ruderman & Simpson (2003) was applied to the non-dispersive 
Alfvenwave, it is shown in the present work that such a transformation is equally applicable 
to a more general, dispersive Alfven modes. 
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In the present work, it is assumed that the charge on dust grains is constant. The 
basic formulation is discussed in Sec. II. It is assumed that the electrons and ions are 
well coupled to the magnetic field and essentially move together. This will restrict the 
applicability of the result to the regions of molecular cloud where neutral number density 
n n < 3 x 10 9 cm~ 3 (Desch & Mouschovias 2001; Tassis & Mouschovias 2005). Assuming 
such equality, the induction equation can conveniently be written in the neutral frame. The 
resulting induction equation is shown to contain the Hall term, due to the relative drift 
between the plasma and the dust grains. It is known that this feature is a generic property 
of the multi-component plasmas (Wardle 1999). The Alfvenpump wave is discussed in Sec. 
Ill and linearized equations are derived. Applying a linear transformation, a set of equation 
with constant coefficient is derived. In Sec. IV the dispersion relation is numerically solved. 
The dependence of the parametric instability on the pump wave amplitude and on the ratio 
of the pump wave to the dust-cyclotron frequency is discussed. In Sec. V, application of 
the results to the molecular cloud is discussed, and conclusion and brief summary of the 
results are presented. 



2. FORMULATION 

We shall assume a weakly ionized molecular cloud consisting of plasma particles, i.e. 
electrons and ions, charged grains and neutral particles. The dynamics of such a partially 
ionized mixture is described by the multi-fluid system. The basic set of equations describing 
the dynamics of such a system is as follows. The continuity equation is 

^ + V-( PjVj ) = 0. (1) 

Here, pj is the mass density and v,- is the velocity of the various plasma components and 
neutrals. We note that the plasma velocities Vj are written in the neutral frame. The 
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momentum equations for electrons, ions, grains and neutrals are 

/ , v,xB\ 

= - qj rij (E + J - Pj u jn Vj , (2) 

Pn jV n = -VP + Pi V 3n Vj • (3) 

Here qj rij (E' + Vj-xB/c) is the Lorentz force and E' = E + v n xB/c is the electric field in 
the neutral frame with E and B as the electric and magnetic fields respectively, rij is the 
number density, and, j stands for electrons (q e = — e), ions (qi = e) and dust qa = Z e with 
Z as the number of charge on the grain and c is the speed of light. 

The collision frequency is 

< cry >j 

Vjn — IjnPn . Pn ■ y*) 

m n + rrij 

Here < av >j is the rate coefficient for the momentum transfer by the collision of the j 
particle with the neutrals. 

The ion-neutral and electron-neutral rate coefficients are (Draine et al. 1983) 

< av > in = 1.9xl(r 9 cm 3 s _1 
< av > en = 4.5xl0~ 9 (J^j 2 em's" 1 . (5) 

Adopting a value of mi = 30 m p for ion mass and m n = 2.33 m v for mean neutral mass 
where m p = 1.67 x 10~ 24 g is the proton mass, the ion neutral collision frequency can be 
written as 

i~n n n n < a v >j n in— 10 -l tr\ 

^in = Pnlin = ; = 1-4 X 10 Tl n S . (6) 

For very small grains (~ 3 — 3000 A), the grain-neutral collision frequency is of the same 
order as the ion-neutral collision frequency, given by equation (6). However, for larger 
grains, the grain-neutral collision rate can vary between 10~ 10 to 10 _5 s _1 for sizes ranging 
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between a few Angstrom to a few microns. This can be seen if we write the collision rate as 
(Nakano & Umebayashi 1986) 

<™> dn =2.8xl0- 5 T 3 |a 2 _ 5 , (7) 

where T 30 is the gas temperature and a_ 5 is the grain radius in the units of 30 K and 
1CT 5 cm respectively. 

We shall define the mass density of the bulk fluid as p rs p n . Then the bulk velocity is 
u pa v n . The continuity equation (summing up equation (1)) for the bulk fluid becomes 

| + V-H = 0. (8) 

The momentum equation can be derived by adding equations (2) and (3) 

du JxB 

As noted above, we shall assume that the electrons and ions are well coupled to the 
magnetic field. Defining 

ft = ^ , (10) 

as the ratio of cyclotron u C j = qjB/rrijC to the collision Uj n frequencies, we shall thus 
assume that fa >> fa >> 1, i.e. plasma particles are tied to the magnetic field and move 
together v e ~ v*. Then making use of the plasma quasi- neutrality condition, n e = rij + Zrid, 
expression for the current density J = e (rij — n e ) v e + Z e n d v d yields 

v e = ^-+v,. (11) 

It is known that the grain size may vary between few microns down to few tens of atoms in 
the interstellar medium. In fact, the presence of very small grains like MRN distributions 
(Mathis et al. 1977), and polycyclic aromatic hydrocarbons (Puget et al. 1985) can play 
an important role in the dynamics. Since small grains are numerous and their charge may 
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vary between ±1 and (Nishi et al. 1991), grains could be coupled to the magnetic field 
either directly, i.e. (3d ^> 1 or indirectly, i.e. even when p d <^ 1 but (rii/rid) (3 d f3i>l (Desch 
& Mouschovias 2001). Following Ciolek & Mouschovias (1993), we shall express grain 
velocity v d in terms of electron velocity as 

v<=ife v " (12) 



where 

= 

Thus, equation (11) becomes 



1 + Vnd 



ft- (13) 
1 - O) T 



Zen d 

Taking curl of the electron momentum equation (2) and making use of Maxwell's 
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equation, in the (3 e 3> 1 limit the induction equation can be written as 

(15) 



dB „ 



(uxB) - ^tS jxB 

Zen d 



We shall note that in the absence of dust, Jx~B/Z en d — > 0, and, Hall term will disappear in 
the induction Eq. (15). The ideal MHD description of the two component plasma assumes 
that the relative drift between electrons and ions are absent, i.e., v e = Vj. The reason for 
the Hall effect in such a plasma is due to the relative drift between the plasma particles, 
since E + v e xB/c = 0, can be written as E + VjXB/c = JxB/en e in the presence of the 
drift. However, in the presence of charged grains, or for that matter in any multi-component 
plasma, since n e ^ rij, even when relative drift between the plasma particles are absent, 
i.e., v e = Vj, owing to the presence of the third, charged component, Hall effect will always 
be present. This can be seen from Eq. (11). The Hall effect disappears in such a plasma 
if in addition to the absence of relative drift, one also demand Z rid <C n e . In molecular 
clouds generally Z rid < n e and thus, Hall drift is always present. Therefore, the Hall effect 
in a multi-component plasma can be caused either by the relative drift between the plasma 
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particles or due to the presence of the charged grains. This is the reason why Hall MHD 
description of a multi-component plasma could be the proper description of the magnetized 
protoplanetary disks (Wardle 1999; Balbus & Terquem 2001). 

The ratio of the Hall to the convective term in the induction equation (15) will shed the 
light on the relative importance of the two terms and the role of the grain magnetization. 
The ratio of the two terms can be written as 

^v^lVxBI-^M^d+e) (16) 
V A A-KZen d p d L u cd 

where u has been replaced by the Alfven speed Va = B / y/A n p and Si = VAi/oj ci is the 



ion-inertial scale. Here L is the characteristic scale length of interest and Va% = B / y/A n pi 
is the ion Alfven speed. Generally Si <C L and Hall term will be insignificant in comparison 
with the convective term. However, in the present case Hall term can become comparable 
to the convective term. We note that even when grains are not well coupled to the magnetic 
field, i.e. <C 1, Hall term may still compare with the convection term. Clearly then, Hall 
effect is important in a weakly ionized medium. 



3. THE PUMP WAVE 

One can investigate the dusty plasma dynamics with the help of Eqs. (8), (9) and (15) 
along with an isothermal equation of state. We shall assume that all physical quantities 
depend on z only in the presence of an uniform background magnetic field B = (0, 0, B). An 
exact solution of Eqs. (8), (9) and (15) are finite amplitude circularly polarized Alfven wave 

B x = A cos (j) , B y = A sin , 
U x = U cos0, U y = U sin0. (17) 
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Here <p = ko z — uj t. The wavenumber k and frequency uo of the pump wave are related 
by the following dispersion relation 



ul = hlVl (1±^-J , (18) 

where ui c d = (pd/p)(l/l + 6) Z eB/m d c is the modified dust-cyclotron frequency. Since 
Pd/p = 0.01 in molecular clouds, irrespective of the grain magnetization, the modified dust 
cyclotron frequency will be much smaller than the unmodified dust cyclotron frequency. 
The amplitudes A and Uq of the pump waves are related by 

^-(tfX-s)" 1 - ™ 

We shall assume now that the steady-state background consists of the unperturbed as well 
as the circularly polarized pump waves, i.e. B = (B x (z), B y (z), B) and u = (U x (z), U y (z), 0) 
with a constant density. The linearization of Eqs. (8), (9) and (15) yields 

d5p 



dt 



+ V - (p5u + 8pu) =0, 



d6u + (u + 5u.V) (u + 5u) = -^ + ^VP + i : ^JxB + Jx5B-^JxB 

Vx 



dt p p 2 cp \ p 

\ x An x H -I n ynH - . „ 

Zen d \ p 



(20) 



dt 

The above set of Eqs. (20) is supplemented with an equation of state 5P = C 2 a bp. In what 
follows, we shall also assume that all the perturbed quantities depend on z only. Then 
writing above set of equations in the component form 

0, 



dSux Su dU x 
dt z dz 
d5u y ^ dU y 
dt z dz 

d5u z C 2 d5p 1 f dSB dSB y \ 5p d(B 2 + B 2 ) 





d5p 


dSu z 




m + 


^ dz 


B ( 


'd8B x 


bp dB x 


An p \ 


v dz 


p dz 


B 


}d5B y 


5p dB y 


4 7T p 


V dz 


p dz 



dt p dz 4 7r p \ dz dz J 8 n p dz 
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dbB x . dB x „ cM-u^ „ <9fo/ 2 / <9 2 55,, 1 <95p 
+ 5u z ^r- = B z — B x — ha' 



dt dz dz dz \ dz 2 p dz dz 

dbB y +Su dBy =B 9Su 1 _ B _ a f ggg? _ I (21) 

dz dz y dz \ dz 2 p dz dz J 

where a = cB (1 + 0) /(An Z en d ). While writing above equations, 5B Z = has been 
assumed. Eqs. (21) have periodic coefficients and belongs to a general class of equations 
known as Floquet equation. The Floquet theorem prescribes the form of the solution 
for the differential equation with periodic coefficients. However, Floquet theorem also 
states that there exists a linear transformation that can reduce the system of ordinary 
differential equations (ODEs) with periodic coefficient to a system of ODEs with constant 
coefficient. For a non- dispersive MHD case, such a transformation has recently been 
proposed (Ruderman & Simpson 2003). In the present case, applying same transformation, 
it is shown that the system of Eqs.(21)can be reduced to ODEs with constant coefficient. 
Defining 

5B + = SB X cos + b By sin , 5B_ = SB X sin — 5B y cos , 

bu + = 5u x cos + 8u y sin , bu_ = 5u x sin — 5u y cos , (22) 

Eqs. (21) can be reduced to following set of equations with constant coefficient. 

dbp dbu z = Q 

dt dz 

d5u + . B fd5B + , ro 

— f- - u bu„ = - — -± + k SB. 

dt An p \ dz 

96u - r , rr S B ( dSB ~ , Sn \ i B A r 

-— — + cj bu+ - k Uq bu z = — k 5B + + k - bp , 

dt An p \ dz J An p l 

dbu z C 2 dbp A dbB + 

dt p dz An p dz 

dbB + rn n f dSu + , s \ fd 2 bB„ n1 dbB + 

-Lu bB_ = B[ -j-^ + k 6u-)-a[ -2k + k 2 bB_ 



dt \ dz J \ dz 2 " dz 

dbu z a k A dbp 
dz p dz ' 

d6B_ ^ „{dbu_ , r \ fd 2 bB + , d6B- l2p „ 
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+k A 5u z , (23) 

One notes that the x and y components of the induction Eq.(Eqs. (21)) are symmetric 
with respect to the Hall terms (with coefficient a) whereas in the Eq. (23) such symmetry 
does not exist. The reason for this loss of symmetry is due to the form of 5B + and 8B_. 
The terms proportional to the density fluctuation in the induction Eq. (21), cancels out for 
5B_. This results in the asymmetric equations for 8B + and 8B_. 



DISPERSION RELATION 



The above set of equations has constant coefficient and thus, one can Fourier analyse 
the spatial and temporal dependence of the fluctuations as ~ exp (iujt — ik z). Denoting 
uj = uj/uj , k = k/ko, F± = 1/(1 ± uj/uj cd ), (3 = C 2 s /Vl and C 2 = k 2 (3F following 8 th order 
dispersion relation is derived from Eq. (23). 



uj + a 7 uj + . . . + a\UJ + a = . 

The coefficients aj, . . . a are given in the appendix. The transition from stability to 
instability will proceed through uj = 0. In the u — > limit, 

a 



(24) 



uj = — - 



ai 



(25) 



In the long wavelength limit retaining only ~ 0{k), 0{k 2 ) terms in the coefficients ai and 
a one may write 

A 2 



kF± 



= -2 



B 2 



1 + 



UJ 



kF+ 



P [1 + F± 



UJ 



-F 4 



A 



2 1 



1-F4 



^0 \ 
Ucd) 



(26) 



(27) 



uj cd I B 2 

Recall that for the left-circularly polarized pump waves, F + = 1/(1 + uj /uj cd ) and hence, 
from Eq. (25), a = 0. Thus in view of Eq. (27), one should anticipate that the instability 
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will disappear in the vicinity of k — 0. As will be shown later, the numerical solution of full 
dispersion relation, Eq. (24) indeed supports this conclusion. 

For the right-circularly polarized pump waves, when F_ = 1/(1 — cuo/^cd), the growth 
rate can be written as 



As is clear from above Eq. (28), the growth rate is inversely proportional to the factor 
(1 — uo/uJcd)- This implies that near u>o ~ u c d, when pump is operating near the modified 
dust-cyclotron frequency, the instability can grow resonantly. Therefore, the growth 
of the parametric instability is quite different for the left and right-circularly polarized 
pump. Whereas for the left-circularly polarized waves, the instability does not exist in the 
neighborhood of k = 0, for the right-circularly polarized pump, the instability may become 
large near k — 0. Clearly, for the right circularly polarized pump, unbounded growth of 
the instability is not possible. The linear approximation in which the dispersion relation 
(24) has been derived breaks down for any dependent physical variable / once Sf < f. 
Therefore, in the vicinity of resonance nonlinear terms needs to be considered. This task is 
beyond the scope of the present work. 

Now the dispersion relation (24) is solved numerically. In Fig. 1(a) the normalized 
growth rate is plotted against the normalized wave number for different values of uj^jujd. 
The left circularly polarized pump is assumed. The instability does not exist in the vicinity 
of k — > and k = 0.2 is the threshold of the instability. It is in conformity with the above 
discussion where in the vicinity of k — 0, wave does not grow. This is caused by the 
dissipation of the long wavelength fluctuations by the plasma-neutral collisions. With the 
increase in uo/uj c d, the growth rate decreases. As is clear from the pump wave dispersion 
relation (18), the increase in uo/u c d implies the increasing importance of the Hall term 
(J x B), which appears due to the relative drift between the plasma and the neutral. This 
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drift is caused entirely by the collisonal momentum exchange. Therefore, the increase 
in u /u c d implies the increased dissipation of the free energy. Hence with the increasing 
L} /u cd , one would expect a decrease in the growth rate. With decreasing u /u cd , the growth 
rate is due to nondispersive Alfvenpump and will correspond to a non-dissipative, ideal 
regime. Hence one sees the saturation of the growth rate with decreasing with decreasing 

^0/ ^cd- 

In Fig. 1(b) the growth rate of the parametric instability is shown against the variation 
of the pump amplitude A /B. Recall that B is the amplitude of the constant magnetic field 
in the z direction and, A is the amplitude of the transverse x — y fields. The growth rate 
is not very sensitive towards smaller wavelength, i.e. after k/k > 0.5. However, the onset 
of instability is slightly delayed when Aq/B — 1 than when Aq/B = 0.1 implying that the 
effect of collision is less pronounced for a stronger B z field. 

In Fig. 2 both real and imaginary frequencies are shown for Aq/B = 1 and ujq/uj C( i = 1. 
We see that from the onset until the demise of the instability, the real part of the frequency 
is constant. The part of free energy that is responsible for setting the normal oscillation in 
the system has been diverted towards the wave growth. If one superposes the imaginary 
part on the real frequency, then one gets an almost constant uo r . This should be expected 
since the total energy available to the system is due to the Alfvenpump, and thus, the 
growth of the fluctuation occurs at the expense of the oscillatory real part. 

The growth rate is not sensitive to the plasma (3. Unlike in the non- dispersive, 
non-dissipative case where the increase of (5 leads to the increase in the growth rate 
(Goldstein 1978), in the present dissipative, dispersive case, plasma (5 appears to have no 
role. In the non-dissipative case the increase in (3 implies the unlikelihood of compressional 
wave excitation, the role of plasma f3 in the dissipative case is unclear. The dispersive 
nature of the pump wave is due to the collisonal effects and thus, f3 is not directly related 
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to the growth of the waves. 

In Fig. 3 (a) the growth rate and in Fig. 3(b) corresponding real part of the frequency 
is shown for the right-circularly polarized pump. When Uo/u> c d = 0.9, the growth rate 
becomes very large. The free energy is resonantly pumped into the fluctuations with 
increasing u Q /u cc i. The physical system behaves like a driven oscillator. The resonant 
driving is indirectly related to the neutral-plasma collisions. The relative drift between the 
plasma and the dust causes a Hall field over the dust-cyclotron time. If the Alfvenwave 
propagation time Uq 1 becomes comparable to the dust-cyclotron time, the energy is freely 
fed to the fluctuation by the pump to the plasma particles. Resulting free energy to the 
fluctuation causes the large growth rate. This behavior has been analytically predicted in 
the limiting case by Eq. (28). It should be noted from the corresponding curve in Fig. 3(b) 
that the real frequency shows a sharp decline for ujq/u C( i = 0.9. This suggest that almost all 
the free pump energy is resonantly used in the fluctuation growth. Similar behavior is also 
noted for uJo/uj cd = 0.8 and u} /u cd = 0.01. Since the instability growth rate in this case is 
smaller than when uoq/uo C( i = 0.9, the part of the available free energy remains in the real 
part of the frequency. This is in agreement with the well known behavior of the oscillators 
near resonance although present system is more complex. 



5. APPLICATION AND SUMMARY 

The grains are important charge carriers in the molecular clouds and their presence can 
significantly alter the dynamics of the low frequency Alfvenwave in the neutral medium. 
Observations of molecular clouds suggest a magnetically threaded supersonically turbulent 
environment, and the generation of magnetohydrodynamic waves in such an environment is 
inevitable. Thus taking the upper value of the pump wave frequency (Folini et al. 2004), 
ujq ~ 10~ 4 yr _1 , the growth rate of the left-circularly polarized wave 0.3 u Q suggests that 
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the parametric instability of the Alfvenwave could be relevant to the onset of turbulence. 
The Alfven wavelength lies in the range of 0.07 — 0.35pc(Folini et al. 2004) and since the 
maximum growth rate occurs at k/k = 0.5 (Fig. 1(a)), the parametric instability of the 
left-circularly polarized mode will operate between 0.1 — 0.7 pc. 

To check whether the resonance condition can prevail in the molecular clouds, we 
shall assume that the grains are magnetized. Since smaller grains are more likely coupled 
to the magnetic field than the large grains, we take m g = 10~ 18 g. For a typical magnetic 
field B ~ 10//G, and taking p^j p = 0.01 and = 10, one gets uj c d ~ 10~ 4 yr _1 which 
is comparable to the pump frequency. Therefore, the right circularly polarized wave can 
resonantly excite this instability. Furthermore, it can operate on very long wavelengths. 

To summarize, parametric instability of a weakly ionized medium may play an 
important role in exciting the turbulence in the interstellar medium. Following is the 
itemized summary of the present work. 

1. The parametric instability of the weakly ionized collisonal medium is studied in the 
present work. It is shown that the collision causes the dispersive nature of the pump waves. 

2. Using a linear transformation, the linearized equations with periodic coefficients 
are reduced to a form with constant coefficient allowing the validity of the normal mode 
analysis at all wavelengths. 

3. The instability is sensitive to the change in the ratio of the pump and the cyclotron 
frequencies and weakly sensitive to the ratio A /B. The growth rate is not sensitive to the 
plasma (5. 

4. The instability can become an order of magnitude larger for the right hand circularly 
polarized pump particularly near the resonance, i.e. when ojq ~ uj C( i- 
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A. Appendix material 
Coefficient of the dispersion relation (24) 
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Fig. 1. — The dependence of the growth rate on the ratio of the left circularly polarized pump 
wave to the dust-cyclotron frequency is shown in Fig. 1(a) for the corresponding physical 
parameters in the box. In Fig. 1(b) the dependence of the growth rate on the amplitude of 
the pump wave Aq/B is shown. 
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Fig. 2. — Both real and imaginary part of uj/uq is plotted against k/ko for the plasma f3, 
ujq/uj C( i and A /B shown in the box. 



-27- 




Fig. 3. — The dependence of the growth rate on the amplitude of the right handed circularly 
polarized pump wave is shown in Fig. 3(a) and corresponding real part of the frequency is 
shown in Fig. 3(b) for different values of uj/uo c d. The value of plasma j3 and Aq/B is 0.1. 



